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Abstract 

We consider a basic stochastic particle system consisting of N identical particles with isotropic 
fc-particle synchronization, k > 2. In the limit when both number of particles N and time 
t — t{N) grow to infinity we study an asymptotic behavior of a coordinate spread of the par- 
ticle system. We describe three time stages of t{N) for which a qualitative behavior of the system 
is completely different. Moreover, we discuss the case when a spread of the initial configuration 
depends on N and increases to infinity as N ^ oo. 

1. Introduction. A wide class of probabilistic models can be interpreted as stochastic particle 
systems with a synchronization-like interaction. In general terms, the matter concerns a special 
class of jump Markov processes x{t) = {xi{t), . . . ,XN{t)) evolving in continuous time and tak- 
ing their values in M^, with generators of the following symbolic form L = Lq + . The vari- 
able Xi G M is interpreted as a coordinate of the particle i. A free dynamics Lq corresponds 
to independent movements of the individual particles, and corresponds to synchronizing jumps 
X = {xi, . . . ,xn) — > x' = {x[, . . . , x'jy) which, according to many preceding papers [2-4, 7], fit to 
the following general rule 

{x[, . . .,x'^} C {xi, . . .,xn} , {x'l, . . . ,x^} / {xi, . . .,xn} . 

First mathematical papers on stochastic synchronization dealt with the case N = 2 (see, for exam- 
ple, [1]). Many-particle systems (A^ > 2), studied till now, have different forms of the synchronizing 
interaction. So [2-4] considered pairwise interactions, while in [7] a three-particle interaction was 
studied. In the present paper we consider a general k-particle symmetrized interaction which will be 
defined in terms of synchronizing maps. 
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Figure 1: Correspondence between maps TTki,.. .,ki and oriented graphs 



When we analyze a collective behavior of such particle system we face with superposition of two 
opposite tendencies: with the course of time the free dynamics increases the spread of the particle 
system while the synchronizing interaction tries to decrease it. The aim of this paper is to study 
a qualitative balance between these two tendencies in the limit N ^ oo, t = t{N) —>■ oo. Our 
main result (Theorem O determines three phases of different behavior of the system depending on the 
growth rate of t{N). Similar result was obtained in [4] for a system with two types of particles and a 
pairwise interaction between different types. 

Main results of this paper were presented by the author at PTAP-2006 (see [8]). 

2. Synchronizing maps. Let us enumerate the particles and let xi, . . . , be their coordinates. 
Let us fix a natural number k >2 and integers ki > 2, . . ., ki > 2 such that A;i + • • • + fc^ = k. We call 
the sequenced collection {ki, . . . ,ki) a signature and will keep it fixed throughout the present paper. 
Let T be a set of all sequenced collections {ii, . . . consisting of k different elements of the index 
set Mn :={!,•••, N}. On the set I we define a map -Kk^^.^^ki '■ {ii, ■ ■ ■ ,ik) ^ (fij ■ ■ ■ , ^i) , where 



lO 



{i2, ■ ■ -Jki), 

(ifc 1+2) • • • ) ^A;i+fc2 j; 



91 = n, 



1 



92 — iki+1, 



^ 2 





It is useful to associate with the map vTfc^ ... an oriented graph as shown on the Figure[T] 
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Consider also a map ^ ■ ^ 



which is defined as follows: 



j Xm, if m ^ (Zi, . . . , 

ym = < / , ^ (1) 

[ , if me (ii,..., Ik), me Tj . 



It is natural to call this map a synchronization of the collection of particles ^ , . . . , Xj^. , corresponding 
to the signature {ki,. . . , ki). It is evident that the configuration Jk^'"''l^''x has at least ki particles 
with coordinates that are equal to , . . . , at least ki particles with coordinates that are equal to Xg^ . 



3. Multidimensional Markov process with synchronization. We define a continuous time homo- 
geneous Markov process x{t) = {xi{t),... , xjv(t)), t > 0, on the state space by means of the 
following formal generator L = Lq + L^. Free dynamics Lq corresponds to independent random 
walks and is chosen as 

JV 

(Lo/)(x) = a^ / (/(x + 2;ei)-/(x))p(dz), zeR, = (0, . . . , 1^, . . . , 0) , 

i=i 

i.e., independently of the other particles each particle i waits an exponential time with parameter 
a > and performs a jump of the form Xi ^ Xi + z, where z is distributed according to the law 
p{dz), which is common for all particles. We assume that the distribution p has a compact support 
and is nontrivial in the following sense: 62 := / x^p{dx) > 0. Denote also a := J xp{dx). 
Synchronizing interaction Lg has the form 

"^^"^ ^ TVW ^ -^s,(ii,...,ifc)/) (2) 
(n,...,ife)ex 

where AtM := N{N - 1) • • • (AT - A; + 1) and 

{x) = f ^x) - fix), x = (xi, . . . , x^v), o) 

for bounded continuous functions /. In other words, independently from the free dynamics at epochs 
of a Poissonian flow with parameter 5 > we choose with an equal probability one element 
{ii, . . . ,ik) from the set X, and synchronize the particle configuration xi, . . . , xjv according to the 
map Jk^'"''lf ■ Hence, the interaction, corresponding to the given choice of ,ik), consists in 

instantaneous and simultaneous transfer of particles of each group Xk^^ \-kj-i+i, ■ ■ ■ ,Xki-\ \-kj to 

the point with the coordinate Xg^ = Xfe^+...+jfc^._j+i, j = 1, 1. 
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Since in the sum ^ all summands have the same weight, the system under consideration belongs 
to the class of models with mean-field interaction. 



4. Main results. Consider a function V : 



pN 



N(N 

It is easy to see that this function coincides with the empiric (sample) variance 5^, where 

N N 
m=l m=l 

Lemma 1 (LM) (x) = aa for all x = {xi, . . . ,xn)- 

Lemma 2 If N ^ oo then we have LV = abo ; r V for some x > 0. 

N{N - 1) 

I 

It follows from the proof that x = A;| — k. 

i=i 

We are interested in the following means: fiN{t) '■= EM{x{t)) and Rwit) := EV{x{t)). In the 
below statements we always assume that ^ oo. 

Theorem 1 For any t > lim ^ ^ 1^n{0) _ Moreover, for any function t{N) oo the 
jollowing convergence holds t{N) ^ " 

Theorem 2 Assume that supjy I-RAf(0)| < oo. There exist three time scales t = t{N), where 
Riy{t{N)) has completely different asymptotic behavior 

• //^^ ^ 0, then RN(t(N)) ~ ab2t(N). 

• Ift{N) = cN^, then RN{t{N)) ~ ab2{6>i)-\l - exp {-6>€c))N^. 

• ^^"t^ ^ ^'^^^ RN{t{N)) ~ ab2{Sx)-^N^. 

In [4] similar consecutive stages were called correspondingly a phase of initial desynchronization, a 
phase of critical slowdown of desynchronization and a phase of final stabilization. We see that the 
first phase does not contribute to the asymptotic behavior. 

Theorem 3 Assume that Riv(O) — > oo as N —>■ oo. Then the additional condition — } is 

^ ' t{N) 

sufficient to ensure the validity of the corresponding statements of Theorem^ 
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5. Proof of LemmaHl A straightforward calculation shows that {Lq M) {x) = aa, hence we need 
only to prove that {L^M) (x) = 0. Let us show that this fact follows from the symmetry of the 
synchronizing interaction. Indeed, if the signature of interaction (fci, . . . , fc/) is given we can fix some 
sets 

Bj G Mn, \Bj\ = kj, [j = Bj, n Bj, = {ji 7^ is) 

and define X^^'"'^^' := ^{ii, . . . ,ik) ■ {gj}L)T° = Bj Vj = 1, z| , where and r° are determined 
by the map vrfci,...,^; applied to (ii, . . . , i^). To finish the proof it is sufficient to show that 

E (^^ ° '^t::^ -m)=o, (4) 

since = S . It is easy to see that 

M - ^(-) = ^ E E - -'^) ■ 

i=i her° 

Substituting the r.h.s. of (O into the sum we see that, since in ^ gj covers the whole set Bj, for 
any pair of indices u,v G Bj we have exactly one difference {xu — Xy), when gj = u, and exactly 
one difference {xy — when gj = v. Hence the sum (01) is equal to zero. 

6. Proof of Lemma m Define a function /m,n : — > M by the formula 

So V{x) = fm,ri{x)- It is Straightforward to check that Lq fm,n{x) = (^N{N — 1)^ 2ab2 and 

hence (Lq V) (x) = 062- 
We have 

L,V = 6 ^ ^s,(n,...,ifc)/m',n' • (6) 

m'<n' («!,.. .,ifc)Gl 

Consider a summand -/^s,(n,...,ifc)/m',n'- By ([B the map jj.^^' " '^^^ transfers a particle with index m' to 
the point having the coordinate Xm, where 

/ m', if m' ^ (ii, ...,ik), 

771 = < (7) 

[ gj, if 777,' G (n,. . .,ik),m'€ Tj 
Consequently, i^s,(n,...,ifc)/m',n' = /m,™ — fm',n' with some 777 and n, which are not necessarily 
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different from m' and n' . Hence, there exist such anm{N) € M that 

{L,V) {x) = 6Y^ amn{N)U,,n{x). (8) 

Our goal is to show that the coefficients amn{N) do not depend on m and n and that 

amn{N) = -^^(^^ (9) 

for some constant x > 0. Let us fix some pair {m,n} and calculate amn{N). When we choose 
a collection (ii, . . . , i^) for the syncfironization, in that way we choose gi, . . . ,gi and r°, . . . , F^. 
Denote G = {gi, . . . ,gi}. If the pair {m, n} is given then the set I can be divided into four non- 

3 

intersecting parts as follows: Z = ^T™'", where 

■w=0 

2-m,n . ^ {(ii, . . . : {m,n} n (UjTj) = 0} , 
Jf'": = { (ii,..., ifc) : {m,n} n(u,r°) /0} , 
2-m,n . ^ {(ii, . . . : {m,n} C G} , 
jrn,n , ^ { (i^ , . . . , ) : | {m, n} R G[ = 1} \T™'" . 



In each sum ^^""^ := ^ ^ (/m',n' o " /m',n') we pick out only sum- 



mands that contain the function fm,n and a coefficient in front of it will be denoted by amn{N)- By 

3 , 

representations Q and ([D we can write anm{N) = X] flrnn(A^)- 

0) If (ii, . . . , ifc) G I™'" then the particles with indices m and n are fixed points of the map 

Uence a^^^(iV) = 0. 

1) If (zi, . . . , ifc) G I™'", then the summand fm,n can be presented in the sum X]*^^^ only with the 
sign " — " and only in the case when (m', n') = (m, n). A total number of collections (zi, . . . , ik), 
that belong to the set If"'", is equal to |2f '"| = iYl*^! - (iV - 2)[*--'l {N - {k - l)f\ Therefore, 

II) For definiteness let us fix some i,j € {1,...,/} and consider a subsum of ^^^^ taken over 
such subset of 2^'" that m = g-i, n = gj. Under the action of the map J^^'''l^^ each of function 
where u G Fj, G Fj, turns into the function Thus in the above subsum we find kikj summands 
fmn with the sign " + " and only one summand (corresponding to the case (m', n') = (m, n) ) with 
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the sign " — ". Hence, 



in) Fix i E {1, . . . , /} and consider in ^^^^ a subsum taken over such {ii, . . . , i^) that m = gi G 
G, n ^ G. A number of collections (ii, . . . , i/^), satisfying to this assumption, is equal to (A^— 2)['^~^1. 
Under this assumption each of ki functions fun, where u G Fj, turns into fmn- Changing the roles of 
m and n, we conclude that 

aE(iV) = 2^(fc.-l) ^""-^f . 
1=1 

As is easy to see the values of a^mn{N), amn{N) and amn{N) do not depend on m, n. After some 
calculations we obtain the following formulae 



2{k - 1){N -k) + {k- l){k + 1-1) 
A^(A^- 1) 

2{k-l){N -k) 



"^"^ ^ N{N-l) 



N{N -1) 
Summing these values together, we get 

a^rtliN) + a^^liN) + a^^liN)-. 
and the statement ^ is proved with x = kj — k > 0. 



N{N-1) 



7. Proofs of theorems. Our method is similar to the approach which was proposed in [4]. It is based 
on an embedded Markov chain Cj\f{n,UL>) := x{Tn,uj), n = 0, 1, . . . . The chain CAf(n) is defined on 
the same probability space (il, T, P), evolves in discrete time and take its values in the set M^. A 
sequence Ti(a;) < T2{lu) < • • • consists of time moments at which particles make jumps, tq{uj) = 0. 
Evidently, {r„+i — t„},^q is a sequence of i.i.d random variables having exponential distribution 
with the mean value 7jv = (aN + 6)^^. In the discrete Markov chain CAf(^) one-step transitions 
have the following form: with probability Q!77v particle i makes a jump Xj — > Xj + z, where z is 



distributed according to the law p{dz), i = l,N,or: with probability 5^^ a synchronization take place, 
namely, a collection of indices (ii, . . . , ij^) is chosen randomly and then the particle configuration x 
is transformed into Jf^^' '^^^- By the law of large numbers t„ ~ •jN'n as n — > oo, therefore, an 
asymptotic behavior of the particle system x{t),t > 0, can be reduced to asymptotic properties of the 
embedded Markov chain CN{n). We are interested in the following sequences: s{n) := EM {x{Tn)) 
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and d{n) := Ey(x(r„)). Straightforward calculations together with Lemmas [T] and |2] show that 

E (M(C7v(n + 1)) I x{t), t<Tn) = M(CivN) + 7iva«, 

E{V{CN{n + l)) \x{t),t<T^) = V{CN{n)) + 7Nl^ab2-dj^^^^—^V{CNin))y 

Taking expected values of the both parts in the above equations we come to recurrent relations 

s(n + l) = s(n) + 77vaa, (10) 
d{n + l) = d{n) ^1 - IN^ _ + lNab2 (11) 

From ([Tol l the statement of Theorem [T] easily follows. By iterating the equation (ITTI ). we get 

n-l 



d{n)=d{0){l-jN6j;^_ 



+ (7^aM E (l - ^^^]^(]^ J . (12) 



Substituting n = j^^t{N) in (fT2l) and letting go to the infinity, we come to the conclusion of 
Theorem |2l 

Theorem [3] can be derived from (fT2l) by a straightforward calculation performed separately for 
each time stage. 

8. Possible generalizations and perspectives. Undoubtedly the results obtained in this paper re- 
main also true, if the synchronization Lj is taken in the form of a symmetrized polynomial interaction 
of any order, i.e., in the case when is a finite sum (taken over index A; > 2) of fc-particle sym- 
metrized interactions. Apparently, the methods of our paper can be adapted also for other classes of 
the free evolutions Lq , such as independent Brownian motions. 

The stochastic particle system considered here plays a basic role for future study of asymptotic 
behavior of general many-component stochastic systems with synchronization. The statement of The- 
orem |2] put forward a hypothesis that, seemingly, the result of the paper [4] on the existence of three 
phases of collective behavior remains true for a wide class of large particle systems with synchroniza- 
tion and, hence, it is not really caused by the specific nature of the model [4]. As it is seen from the 
present paper, key elements in the future proofs of such results should be some analogues of Lemma|2l 
On this way we expect to have difficulties with anisotropic synchronizations which are interesting in 
a number of important applications. Let us remark that in papers [2,7], where some examples of 
anisotropic interactions were considered, the behavior of particle systems was studied only on a so 
called hydrodynamic scale while in the present paper we consider all possible time scales t{N). 

Unfortunately, an important class of cascade synchronizations [5, 6] can not be considered in the 
framework of the present approach. 
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